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О задаче Коши для уравнения в частных производных первого порядка 
и её приложениях в теории обратных задач 


А. О. Ватульян, Л. С. Гукасян 
(Донской государственный технический университет) 


Исследуется обратная коэффициентная задача для оператора второго порядка в односвязной области с ку- 
сочно-гладкой границей, возникающая в теории колебаний деформируемых систем. Предложен метод реше- 
ния обратной коэффициентной задачи на основе исследования задачи Коши для дифференциального урав- 
нения первого порядка с переменными коэффициентами. Решены прямая и обратная задачи на основе мето- 
да разностных аппроксимации. Приведены результаты реконструкции переменного модуля сдвига различных 
типов, полученные как при точных, так и при зашумленных входных данных. 

Ключевые слова: задача Коши, обратная коэффициентная задача, разностные схемы. 


Введение. Коэффициентные задачи в естествознании, в математической физике, в механике 
деформируемого твёрдого тела — интенсивно развивающийся раздел вычислительной и экспери- 
ментальной механики, требующий основательной теоретической базы, составляет один из важ- 
нейших классов обратных задач. Для этого класса задач можно выделить модели, в которых 
идентифицируемые дифференциальные операторы имеют постоянные коэффициенты и модели, в 
которых требуется идентифицировать неоднородные свойства [1, 2]. 

Заметим, что можно рассматривать две кардинально различные постановки задачи иден- 
тификации неоднородностей в зависимости от способа измерения полей. В первой постановке 
реконструкция искомых характеристик осуществляется по измерению внутренних полей, а во вто- 
рой — по измерению граничных полей в некотором частотном диапазоне. В первом случае обрат- 
ная задача линейна, во втором — существенно нелинейна. Переменность коэффициентов диффе- 
ренциальных операторов не позволяет построить в явном виде общие представления решений 
для соответствующих операторов, как для операторов с постоянными коэффициентами. Если ко- 
эффициенты дифференциальных операторов меняются произвольным образом, то методы реше- 
ния прямых задач опираются либо на аппарат интегральных уравнений Фредгольма второго рода 
(для стержней и пластин) [3] и численную процедуру обращения соответствующих конечномер- 
ных операторов, либо на прямое использование конечноэлементных технологий. В то же время 
анализ первой линейной постановки позволяет глубже разобраться с причинами сильной некор- 
ректности, возникающей при анализе коэффициентных обратных задач [4]. С математической 
точки зрения эта задача сводится к решению некоторой задачи Коши для системы дифференци- 
альных уравнений в частных производных первого порядка [5]. В данном случае известны раз- 
личные способы анализа. Для обыкновенных дифференциальных операторов такие постановки 
рассмотрены ранее — в [6, 7] проанализированы подобные обратные задачи для уравнении пер- 
вого и второго порядков. 

В приложениях часто встречается уравнение в частных производных второго порядка при 
постоянных характеристиках. В настоящей работе рассмотрена задача Коши для такого уравне- 
ния, известного как уравнение Гельмгольца. 

1. Постановка прямой задачи. Рассмотрим установившиеся крутильные колебания стержня с 
переменным модулем сдвига с поперечным сечением $. Уравнение установившихся колебаний с 
частотой имеет вид [8]: 
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(ни. м + роз?и = 0, 
и (1) 
к” вх,’ 
где и — переменный модуль сдвига, р — плотность. 
Будем считать, что поперечное сечение стержня $ есть односвязная область с кусочно- 
гладкой границей /=/ (1. Часть границы / защемлена, на /› действует нагрузка. В этом слу- 


чае краевые условия имеют вид: 


[1 





Й = 0, ы— =Р. (2) 


Приведём краевую задачу (1), (2) к безразмерному виду, введя безразмерные параметры, функ- 
ции и обозначая м, = тах р, р, = тахр, хе 5, 


2_2 
Н 2 _ Р® а 





д=-—,к ‚а = Чат5$ = тахр(А,В), (3) 
Но Мо А,Вс5 
(9%, ) +ки=0. (4) 
и 
и’ =0, —| =р. 5 
/ 9 дп ь Р% ( ) 


Задача (4)—(5) при известной положительной функции 9(х,,х,) представляет собой краевую 
задачу об отыскании и(х,,х,), для которой несложно установить существование обобщённого 
решения при минимальных требованиях относительно гладкости 9 (х,,х,). К сожалению, реше- 


ние такой задачи даже в простой области типа прямоугольника может быть построено лишь чис- 
ленно, на основе либо конечноэлементного подхода, либо разностных аппроксимаций. 
2. Постановка обратной задачи. Рассмотрим задачу, обратную по отношению к сформулиро- 
ванной выше (4)—(5). Поставим задачу об определении функции 9(х,,х,) в области $ по из- 
вестной (измеренной) функции и(х,,х,). Тогда в предположениях о некоторой гладкости (кото- 
рые будут описаны ниже) имеем следующую задачу Коши: 

Чи. +9, + ддичкти =0 


9 а р й 
дп, ^ 
Далее будем требовать выполнения следующих условий: 
И, (7) 
ОП |, 


2) и, и.,Аи — непрерывно дифференцируемы в $, 

3) функции и‚,и, — коэффициенты при производных искомой функции в (6) — одновременно в 
ноль не обращаются, т. е. и’ чи,’ *0. 

Пусть уравнение кривой /, имеет следующую параметризацию: х, =Ф, (5), х, = (5), 
5; Е . Тогда перепишем граничное условие для задачи Коши (6), которое приобретает вид 


и 


9 (4. (5), (5)}5# | 


= р» (5). (8) 
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Введём понятие характеристик системы [9]: 


АХ, _ в ди 
@ "5х, 
(9) 
аХ, ди 
р —= (/ 5  _ 
АЕ ” ОХ, 
Пусть х, (Е), х, (Е) — её решение, которое называется характеристическим. 
Очевидно, что вдоль характеристик выполняется соотношение 
а9 
и. +924» = рт (10) 


Действительно, если д =9д [ЖЖ то по правилу дифференцирования сложной функции имеем 


49 4, д, -а = ди, +9и,. (11) 


= (в НЕВЫ. 
ЧЕ ^” а 
Имеют место следующие утверждения. 
Утверждение 1. Пусть кривая /, не касается характеристик. Тогда задача Коши (6) одно- 
значно разрешима в некоторой окрестности /). 


Проведём из У точки М, е /› характеристику, т. е. решим семейство задач Коши для системы (9): 


ах 
Е =и: Хо =п(5) 


ах 
а | =В (5), 5Е | $1,5, |. 


Обозначим х, (#,5),х, (#,5) — решение задачи Коши (12). В силу (11) исходная задача (6) пре- 


(12) 


образуется к виду. 


9 + ли +к?и =0 

(13) 
в = 2, ($). 
ны 


Таким образом получим, что (13) представляет собой задачу Коши для обыкновенного диффе- 
ренциального уравнения первого порядка. Нетрудно доказать, что 9=ф(Ё,5) — непрерывно 


дифференцируема. Также можно выразить (2,5) через х,,х, в предположении, что якобиан от- 
личен от нуля: 








Ох, ОХ, 
г Г 
5 = 0$ 05 = 1 2 — 0 р 
ОХ, ОХ, и ре 
ОЁ ОЁ |0 


что выполняется, поскольку /› не является характеристикой. 
Утверждение 2. Решение задачи (6) единственно. 
Доказательство осуществляется аналогично [9]. В предположении наличия двух решений 
9,,9› составим задачу Коши относительно их разности д = д, - 9, 
а9 
%& 
которая в рамках условий (7) имеет только нулевое решение. 


+ 9Аи =0, 9 =0, 
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3. Построение решения на основе разностной трактовки. В настоящем пункте рассмотрим 
далее примеры конкретных расчётов для различных областей. Основное внимание уделим зада- 
чам для прямоугольника и круга, причём основные цели исследования таковы: 


1) разработка метода построения функции и(х,,х,) при различных значениях параметра 
к при произвольных законах изменения 9(х,,х,) (монотонных, немонотонных, кусочно- 


непрерывных); 
2) разработка метода построения решения обратной задачи, которая состоит в нахожде- 
нии функции 9(х,,х,) по информации об узловых значениях функции и (х,,х, } при некотором 


значении к, ; 


3) проведение серии вычислительных экспериментов, позволяющих оценить эффектив- 
ность предлагаемых методов. 

Пример 3.1. Рассмотрим задачу решения дифференциального уравнения (6) с переменны- 
ми коэффициентами в области $ =| 0,1 |х| -с,с | со следующими граничными условиями: 


ОИ 
= 


д.94. 
ОХ 


= Х,), 
Ро(х,) — 


=0, (и О: (14) 





Эта задача представляет собой классическую краевую смешанную задачу для неоднородного эл- 
липтического оператора типа Гельмгольца с переменными коэффициентами, и её решение в об- 
щей ситуации (при произвольных законах неоднородности) возможно построить лишь численно. 
3.1.1. Исследование прямой задачи. Прямая задача состоит в нахождении функции 
и (х, ‚Хх, ), удовлетворяющей краевым условиям (14) при заданном законе изменения 9(х,,х,) в 


некотором диапазоне изменения к. 

Один из методов решения основан на методе разностных аппроксимаций. В методе раз- 
ностных аппроксимаций применён пятиточечный шаблон для вторых производных и двухточеч- 
ный — для первых производных (левая разностная производная) [10]. Введём прямоугольную 
сетку с шагами Л, А, по координатным осям х,,х, соответственно. 


Введём обозначение для узловых  переменнын_х® и, =и(Т,-<+Л) и 


9; = 9(1,,-с + Л, ) . Положим шаг аппроксимации по обеим осям одинаковым А, =П, =Р: 


[Ян И 9; ) (мн г И; ; | + 9, ; бы — Чи, , р И р И лв а Ив т [бе > 9, ; ) (“= — И; ; | + В°кти, , =0. 
После использования разностных аппроксимаций получим следующие граничные условия: 
Ин, = Ро;П / 7 + Им-ь, 
И с =И Им =М; ми, Ш, =0 


Г,—С Г,-С+В! 


(15) 


Далее составим систему линейных алгебраических уравнений. Учитывая граничные условия (15), 


которые позволяют исключить граничные узловые точки, получим 
№М-1 М-1 


2..2 
ет О . та а 2 +И, у, + 9+ [Ме РН +7 ки, =0. (16) 


1=1 = 
Система (16) имеет решение, причём единственное. Решая систему (16), находим узловые значе- 
ния функции и, ,. 
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Проверим точность предложенной схемы на примере, когда задача имеет точное решение. 





2 
Например, пусть нагрузка имеет вид р, (х,) =1+2060$ [7 + 3 = | ‚ 9 =1. Точное реше- 


ние задачи строится методом разделения переменных и имеет следующий ВИД 


п (КХ 
иг (хх, ) = В 





+ 0,068 $1 (2,978х, ) со$ [7 + 0,002 $1 (6,203х, ) со$ [7 | | 
ксо$ (к) С С 
Проведена серия вычислительных экспериментов, которые свидетельствуют об эффективности 
предложенной разностной схемы. 

Для иллюстрации на рисунке 1 представлены результаты решения прямой задачи, причём 
сплошной линией обозначено точное решение, а символом «звёздочка» решение, полученное 
методом разностных аппроксимаций, где с =1, В=0,025, М=40, М =80, к=1. 


25 


к 


и (< 1.0) 


01 02 03 04 04 06 00$ 09 
х]| 


Рис. 1. Результаты решения прямой задачи 


3.1.2. Исследование обратной задачи. Целью обратной задачи является восстановление 
неизвестной функции 9(х,,х,) по известным узловым значениям функции и(х,,х, ), получен- 


ным в результате решения прямой задачи методом разностных аппроксимаций, описанным в 
3.1.1. Основной проблемой на этом пути является некорректная задача вычисления производной 
от функции, заданной таблично. 

Для оценки точности разностного подхода введём в рассмотрение относительную 
погрешность: 


Чт; — Эль; 
\ = ‚трах РИ .100%, 
И Эт! 


где 9,,, — точное решение, 9„,, — полученное. 


Рассмотрим случай, когда и=х,, р, =2+х,*. Точное решение при к=1 имеет следую- 
ЩИЙ ВИД: 
2 
Хх 5 
9(х, ‚Хх, ) а и 


На рисунке 2 представлены результаты вычислительных экспериментов полученных при восста- 
новлении различных функций 9(х,,х,). Здесь и далее прерывистой линией будет обозначена 
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восстанавливаемая функция, а точками — узловые значения, полученные при решении обратной 
задачи, Л, = А, = 0,04, с =1, к=1, №М=25, М =50. 





25 1. 
2 1. 
2 ыы 
Е(1.0) ^ 51.0) | 
2. 1. 
2 1. 
=: 1.2 
1. 
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.5 0.59 
х1 
а) 6) 
9(х,х,) =2+х, +х,',У=2.10°% 9(хьх,)=2-х, -х,’,у=1.10°% 
РА. 
= (10) 


В Ва Ва В Ба а Е № 





0.1 02 03 04 053 06 ОГ ОЕ 09 0.1 0.2 03 0.4 05 0.6 0.7 0.8 059 


х[ х1 
в) г) 
9(х,‚х,)=2+яп (пх, ) со$ (пх, ),\ = 0,19% 9(хих,)=2-х, +х,,Уу=3.10° % 


Рис. 2. Результаты решения обратной задачи 


Для аналогии с реальными экспериментами функция и (х,,х, ), используемая как входная 
информация для восстановления 9(х,,х, ), была зашумлена. Вместо узловых значений и(х,,х, ) 
вводились узловые значения и, (х,,х,), такие, что пах и (хи, х, ) — из (х,,х, | =5. Зашумление 


и(х,‚х›) моделировалось с помощью генератора случайных чисел путём добавления к исходным 


узловым значениям малой случайной функции с равномерным законом распределения. 
На рисунке 3 представлены результаты восстановления функции 9(х,,х,) с зашумлен- 


ными входным данными д =4.103,к=1. 

При увеличении величины 6 наблюдается сильный рост погрешности определяемой 
функции, что связано с некорректностью изучаемой задачи. Для улучшения качества реконструк- 
ции необходимо использовать регуляризирующие процедуры, например сплайн-аппроксимации 
при вычислении производных [11]. На основе описанного метода проведены эксперименты по 
реконструкции различных видов неоднородностей: гладкой монотонно возрастающей функции, 
гладкой монотонно убывающей функции, гладкой немонотонной функции. 
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=. 


(10) 


„а 


5% 1.0) 


Вы Бы Ва В Ва Ва Ба Ба 
: ‚ал . к. , 


РНИИ | 





010203 04 03 06 07 08 09 
х1 
6) 
9(хих,)=2-х, -х,’,\ = 6,72% 








ых в 
|. ее 
13 “ Т» 
1.7 и 
1. 
51.0) 1.5 
1.4 
13 
1. 
0.4 02 03 04 05 06 07 05 09 01 02 03 04 05 06 07 08 09 
х1 21 
в) Г) 
9 (х,‚х,) =2+яп (пх, )со$ (пх, ),\ = 8,03% 9(хих,)=2-х, +х,7,\ = 6,63% 


Рис. 3. Результаты решения обратной задачи с зашумленными входными данными 


Пример 3.2. Рассмотрим задачу решения дифференциального уравнения (6), где область $ 
является кольцом. В полярной системе координат, считая, что д = 9(г), и вводя безразмерные 


параметры и функции 





ро? Ь* т 
к =, 0 =9(1),т БОН. 
| (т 9% 
а 9(605 
РЕ ЗЕЕ, = 9 = ) 
Ь 9 
получаем уравнение 
ори |+ки-0 | (17) 


где граничные условия принимают следующий вид: 
и(&)=0, и’(1)=ть. (18) 
3.2.1. Прямая задача. Аналогично предыдущему используем разностные аппроксимации. 
Применим трехточечный шаблон [10] и введём обозначения и, =и(& +1) и 9, =9(& +71). По- 


лучим следующее разностное уравнение (шаг аппроксимации для простоты записи положим 
В =1, Е=&+М): 
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ПВ 22 
(9, в 91 и; е и 1)+ о 2 +И 1+ Е —- и) + Р^КАЫ, =0, 
граничные условия примут следующий вид: 
П=О, Ире ТО ь 
Составим систему разностных уравнений с учётом граничных узловых точек: 


р: |. —(, ти -и.)|- О (м, -и 1) + бк, =0. 
7= 


Для оценки точности предлагаемого метода сравним решение с точным. Так, положим 
2 
9(&)=2+3=”. 
Отметим, что в рассматриваемом случае краевая задача (17) имеет точное решение, выражающе- 
еся через присоединённые функции Лежандра [12]. 


и (5) = 20,65 едепагеР |: _ и ЗЕ? + 3 _1,6Иедепагео |1 _ и ‚ ЗЕ? + } 
На рисунке 4 изображены результаты решения прямой задачи для а = 2 + 3Е*. Здесь прерывистой 


линией изображено точное решение, кружком — результаты, полученные методом разностных 
аппроксимаций (шаг аппроксимации Л = 0,001) при к=1. 





0.10 0.12 0.14 0.156 0.18 0.20 


Рис. 4. Точное и полученное решение прямой задачи 


3.2.2. Обратная задача. В решении обратной задачи также используем метод разностных 
аппроксимаций, описанный в 3.1.1, где в качестве функции и взяты узловые значения, получен- 


ные при решении прямой задачи 3.2.1. 
На рисунке 5 представлен график восстановления функции д =2- 3? , а также результа- 
ты решения обратной задачи по восстановлению функции с зашумлением б =1.10-? (зашумление 


моделировалось аналогично 3.1.2), где р =1.10>. Здесь сплошной линией обозначено точное 
решение, кружком — результаты решения обратной задачи. 
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а) 6) 
У =1.103% У =14,2% 


Рис. 5. Результаты обратной задачи, полученные при входных точных (а) и зашумленных (6) данных 


Выводы. На основе предложенного метода проведена реконструкция переменного модуля сдвига 
различных типов. Полученные данные позволяют утверждать, что при точной входной информа- 
ции такой подход, основанный на методе разностных аппроксимации, эффективен. Об этом сви- 
детельствует величина относительной погрешности вычислений. При зашумленных входных дан- 
ных можно сказать следующее. Предложенная схема работает достаточно эффективно при малых 
степенях зашумления, но, к сожалению, при увеличении зашумления наблюдается значительный 
рост погрешности определяемой функции, что связано с некорректностью изучаемой задачи. 
Необходимо использовать регуляризирующие процедуры. 
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ОМ САЧУСНУ РКОВЕЕМ РОК Е1В$Т-ОКОЕВ РАВТТАЕ ОТЕЕЕВЕМТТАЕ ЕОФЦЧАТТОМ 
АМО ТТ$ АРРИТСАТТОМ$ 1М ТМУЕВ$ТОМ ТНЕОКУ 


А. О. Ма щуапт, (. $. СиКазуап 
(Роп Чае Тесптса! Упмег®Ку) 


Тре тпуегзе соесепЕ ргоМет Юг {пе 5есопа-огаег орегаюг т пе тр! соппесей аотат и те р/есеиле- 
5тоот Боипоаагу ап$та т те {Веогу оЁ Те аеюгта Ме зу%ет игаНоп$ [5 пуезйдееа. Тйе оийоп тепода Юг 
{Пе пуегсе соеИсепЕ рго Мет Базеа оп 5{шаута Саисйу рго ет Гог те Пг5Е-огаег АТегепвйа! едиайоп и уа!- 
аб/е соегГИсгепеЕс /5 опегеа. Во? а!гесЕ апа туег5е рго ет аге 5о№е4 оп те дгоипа о! {пе а!Тегепсе арргохта- 
Воп5 тетод. Тре гесопягисвоп о! те уапа Ме $Веаг тоди/и$ о! уапои$ урез, офатеа аЁ Бо ассигае апа 
по/[5у приЕ дав, 15 гези[ед. 

Кеуигога$: Саиспу рго Мет, пуег5е сое с!епЕ ргоМет, а!Тегепсе спетее. 
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